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ABSTRACT 

We performed Spitzer Infrared Spectrograph mapping observations covering 
nearly the entire extent of the Cassiopeia A supernova remnant (SNR), producing 
mid-infrared (5.5-35 /im) spectra every 5"-10". Gas hues of Ar, Ne, O, Si, S and 
Fe, and dust continua were strong for most positions. We identify three distinct 
ejecta dust populations based on their continuum shapes. The dominant dust 
continuum shape exhibits a strong peak at 21 fim. A line- free map of 21 /xm- 
peak dust made from the 19-23 fim range closely resembles the [Ar II], [O IV], 
and [Ne II] ejecta-line maps implying that dust is freshly formed in the ejecta. 
Spectral fitting implies the presence of Si02, Mg protosilicates, and FeO grains 
in these regions. The second dust type exhibits a rising continuum up to 21 
fim and then flattens thereafter. This "weak 21 /im" dust is likely composed of 
AI2O3 and C grains. The third dust continuum shape is featureless with a gently 
rising spectrum and is likely composed of MgSiOs and either AI2O3 or Fe grains. 
Using the least massive composition for each of the three dust classes yields a 
total mass of 0.020 Mq. Using the most-massive composition yields a total mass 
of 0.054 Mq. The primary uncertainty in the total dust mass stems from the 
selection of the dust composition necessary for fitting the featureless dust as well 
as 70 /xm flux. The freshly formed dust mass derived from Gas A is sufficient 
from SNe to explain the lower limit on the dust masses in high redshift galaxies. 
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Introduction 



The recent discovery of huge quantities of dust (10^ — 10^ Mq) in very high-redshifted 
galaxies and quasars (Isaak et al. 2002; Bertoldi et al. 2003) suggests that dust was produced 
efficiently in the first generation of supernovae (SNe). Theoretical studies (Kozasa et al. 
1991; Todini & Ferrara 2001, hereafter TF; Nozawa etal. 2003, N03) predicted the formation 
of a significant quantity of dust (~ 0.1 — 1.0 Mq) in the ejecta of Type II SNe, and the 
predicted dust mass is believed to be sufficient to account for the quantity of dust observed 
at high redshifts (Maiolino et al. 2006; Meikle et al. 2007). Recently, a model of dust 
evolution in high-redshift galaxies (Dwek et al. 2007) indicates that at least 1 Mq of dust 
per SN is necessary for reproducing the observed dust mass in one hyperluminous quasar 
at z = 6.4. Observationally, the presence of freshly formed dust has been confirmed in a 
few core-collapsed SNe, such as SN1987A, which clearly have showed several signs of dust 
formation in the ejecta (see McCray 19 93 for details). The highest dust mass obtained 
so far for SN 1987A is 7.5 x 10"^ Mq jErcolano et all l2007h . Recent Spitzer and HST 
observations (Sugerman et al. 2006) showed that up to 0.02 Mq of dust formed in the ejecta 
of SN2003gd with the progenitor mass of 6-12 Mq, and the authors concluded that SNe 
are major dust factories. However, from the detailed analysis of the late-time mid-infrared 
observations, Meikle et al. (2007) found that the mass of freshly formed dust in the same 
SN is only 4 x 10~^Mq, and failed to confirm the presence 0.02 Mq dust in the ejecta. The 
aforementioned results show that the derived dust mass is model-dependent, and that the 
amount of dust that really condenses in the ejecta of core-collapsed SNe is unknown. 

Cassiopeia A (Cas A) is the only Galactic supernova remnant (SNR) that exhibits clear 
evidence of dust formed in ejecta (Lagage et al. 1996; Arendt et al. 1999, hereafter ADM). 
The amount of dust that forms in the ejecta of young SNR is still controversial. Previous 
observations inferred only < 3 x 10^'^ Mq of dust at temperatures between 90 and 350 K 
(ADM; Douvion et al. 2001, hereafter DOl) ; this estimate is 2 to 3 orders of magnitude too 
little to explain the dust observed in the early Universe. Recent submillimete r observations 
of Cas A and Kepler with SCUBA founne et all I2OO3I : iMorgan et aD l2003h revealed the 
presence of large amounts of cold dust (~ 0.3 — 2M0 at 15-20 K) missed by previous 
IRAS/ISO observations. On the other hand, highly elongated conductive needles with mass 
of only 10~^ to 10"^ Mq could also explain a high sub-mm fiux of Cas A, when including grain 
destruction by sput tering ( Dwek 2004), though the physicality of such needles is doubtful 
( IGomez et al.ll2005l ). While iKrause et al.l (120051 ) showed that much of the 160/im emission 
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observed with Multiband Imaging Phot ometer for Spitzer (MIP S) is foreground material, 
suggesting there is no cold dust in Cas A. IWilson fc Batrla I (120051 ) used CO emission towards 
the remnant to show that up to about a solar mass of dust could still be associated with the 
ejecta, not with the foreground material. These controversial scenarios of dust mass highlight 
the importance of correctly identifying the features and masses of dust freshly formed in Cas 
A. The Galactic young SNR Cas A allows us to study in detail the distribution and the 
compositions of the dust relative to the ejecta and forward shock with infrared Spectrograph 
onboard the Spitzer Space Telescope. 

Cas A is one of the youngest Galactic SNRs with an age of 335 yr attributed to a SN 
explosion in AD 1671. T he progenit or of Cas A is believed to be a Wolf-Rayet star with 



high nitrogen abu ndance (lFesenll200lh and to have a mass of 15-25 Mq ( iKifonidis et al.ll2000 



Young et al.ll2006l ) or 29-30 Mq ( iPerez-Rendon et al.ll2002l ). The predicted dust mass formed 
in SNe depends on the progenitor mass; for a progenitor mass of 15 to 30 Mq, the predicted 
dust mass is from 0.3 to 1.1 M0 (N03) and from 0.08 to 1.0 Mq (TF), respectively. In this 
paper, we present Spitzer Infrared Spectrograph (IRS) mapping observations of Cas A, and 
identify three distinct classes of dust associated with the ejecta and discuss dust formation 
and composition with an estimate of the total mass of freshly formed dust. 



2. The IRS Spectra and Dust Maps 

We performed Spitzer IRS mapping observations covering nearly the entire extent of 
Cas A on 2005 January 13 with a total exposure time of 11.3 hr. The Short Low (SL: 5-15 
/im) and Long Low (LL: 15-40 ^m) IRS mapping involved ~16x360 and 4 x 91 pointings, 
producing spectra every 5" and 10", respectively. The spectra were processed with the S12 
version of the IRS pipeline using the CUBISM package (Kennicutt et al. 2003; Smith et 
al. 2007), whereby backgrounds were subtracted and an extended emission correction was 
applied. The spectral resolving power of the IRS SL and LL modules ranges from 62 to 124. 

The IRS spectra of Cas A show bright ejecta emission lines from Ar, Ne, S, Si, O, and 
Fe and various continuum shapes as indicated by the representative spectra in Figure 1. 
The most common continuum shape exhibits a large bump peaking at 21 /im as shown by 
spectrum "a" in Figure [H This "21 /xm-peak" dust is often accompanied by the silicate 
emission feature at 9.8 /im which corresponds to the stretching mode. A second class of 
continuum shapes exhibits a rather sharp rise up to 21 /im and then stays flat thereafter. This 
"weak-21 /im dust" is often associated with relatively strong Ne lines (in comparison with 
Ar lines) and is indicated by spectrum "b" in Figure [H The third type of dust continuum is 
characterized by a smooth and featureless, gently rising spectrum with strong [O IV]-|-[Fe II] 
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and [Si II] emission lines as shown by spectra "c" and "d" in Figure [H The spectrum "d" 
shows double line structures that may be due to doppler-resolved lines of [O IV] at 26 /im 
and [Si II] at 35 /xm. Note that the "featureless" dust (spectrum "d" in Fig. [T]) is a class of 
dust, separate from the interstellar /circumstellar dust (spectrum "e" in Fig. [T]) heated by the 
forward shock. The interstellar /circumstellar dust spectrum in Cas A has no associated gas 
line emission. The "broad" continuum (see Figure 7b of Ennis et al. 2006) is a combination 
of the spectra "c" and "e" . The spectrum "c" has contamination from the shock heated dust 
in projection, and for simplicity it is excluded in estimating the masses of the freshly formed 
dust (see §5). The "featureless" dust lacks the gentle peak around 26 fim and also lacks the 
interstellar silicate-emission feature between 9 /im and 11 fim observed in the spectra from 
the forward shock region. Most importantly, the "featureless" dust accompanies relatively 
strong Si and S ejecta lines and mostly from the interior of the remnant (blue region in Fig. 

m). 

We generated a map of the 21 /im-peak dust by summing the emission over 19-23 /xm 
after subtracting a baseline between 18-19 fim and 23-24 /im. The line-free dust map (Fig. 
[2^) resembles the [Ar II] and [O IV] + [Fe II] ejecta-line maps, as shown in Figures and 
and we also find that the [Ne II] map is very similar to the [Ar II] map. The [Ar II] map shows 
a remarkable similarity to the 21 /xm-peak dust map (Fig. [2^ and[2)D), thereby confirming 
this dust is freshly formed in the ejecta. Maps of [Si II] (Fig. [2]1) and [O IV] + [Fe II] (Fig. 
[2t) shows significant emission at the center revealing ejecta that have not yet been overrun 
by the reverse shock (unshocked ejecta). There is also [Si II] and [O IV] + [Fe II] emission at 
the bright ring indicating that some of the Si and 0+Fe ejecta have recently encountered 
the reverse shock. While the bright 0+Fe emission outlines the same bright ring structure 
as the [Ar II] and 21 /im-peak dust maps, the bright part of the Si shell shows a different 
morphology from the other ejecta maps. 

We can characterize the spectra of our three dust classes by using the flux ratios between 
17 /im and 21 /im and between 21 /im and 24 /im. Although the spectra in Cas A show 
continuous changes in continuum shape from strong 21 /im peak to weak 21 /im peak and 
to featureless, we can locate regions where each of the three classes dominates. Figure [2j: 
shows the spatial distribution of our three dust classes where red, green, and blue indicate 
21 /im-peak dust, weak-21 /im dust, and featureless dust, respectively. The flux ratios used 
to identify the three dust classes are as follows where Ix is the flux density in the extracted 
spectrum at wavelength A (/xm): 

1) 21 /im-peak dust: we use the ratio hi/hi > 1+0^21/245 where cr2i/24 is the dispersion in 
I21/I2A over the remnant, which is equivalent to /21//17 > 3.4. The regions with 21 /im-peak 
dust coincide with the brightest ejecta. 
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2) Weak-21 fim dust: we use the ratio 1 — 0-21/24 < -^2i/-^24 < 1 + 0-21/24, which is 
equivalent to 2.3 < /21//17 < 3.4. The regions showing the weak-21 fim continuum shape 
mostly coincide with faint ejecta emission, but not always. 

3) Featureless dust map: we use the ratio /21//24 < 1 — 0-21/24, which is equivalent to 
I21/ In < 2.3. This ratio also picks out circumstellar dust heated by the forward shock, so we 
used several methods to exclude and mitigate contamination from circumstellar dust emis- 
sion. First, using X- ray and radio maps, we excluded the forward shock regions at the edge 



of the radio plateau (iGotthelf et al.ll200ll ). Second, there are highly structured "continuum- 



dominated" X-ray filaments across the face of the remnant which are si milar to the ex terior 



forward shock filaments and may be projected forward shock emission (lDeLaneyll2004l ). For 
our analysis, we excluded regions where there were infrared counterparts to the projected 
forward shock filaments. Third, for simplicity we excluded regions with gently rising spectra 
identified by curve "c" (the spectra which continues to rise to longer wavelengths) in Figure 
[U This type of spectrum is mainly found on the eastern sid e of Gas A w here there is an Ha 



region, the northeast jet, and other exterior optical ejecta (jFesenll200ll ) making it difficult 
to determine if the continuum emission is due to ejecta dust or circumstellar dust. However, 
note that some portion of the continuum in the spectra, "c", is freshly formed dust. We 
finally e xcluded regions where there was a noticeable correlation to optical quasi-stationary 



flocculi (Ivan den Bergh Ill97ll ) which are dense circumstellar knots from the progenitor wind. 



The featureless dust emission appears primarily across the center of the remnant, as 
shown in Figure[21i (blue). The featureless dust is accompanied by relatively strong [Si II] and 
[S III] and [O IV+Fe II] lines, as shown by the spectrum "d" of Figure [H The [O IV] + [Fe II] 
line map (Fig. ^) shows significant emission at the center as well as at the bright ring of 
the reverse shocked material. The [Si II] line map shows different morphology than other 
line maps and the 21 /xm-peak dust map; depicting center-filled emission with a partial shell, 
as shown in Figure [2J This poses the following important question: why is the Si map more 
center-filled than the Ar map? The answer is unclear because Si and Ar are both expected 
at similar depths in the nucleosynthetic layer (e.g. Woosley, Heger & Weaver 2002). The 
relatively faint infrared emission of Si and S at the reverse shock may imply relatively less 
Si and S in the reverse shock. We suspect it is because the Si and S have condensed to solid 
form such as Mg protosilicate, MgSiOs, Mg2Si04 and FeS. In contrast, Ar remains always 
in the gas and does not condense to dust, so it should be infrared or X-ray emitting gas. 
Alternate explanation is that the ionization in the interior is due to photoionization from the 
X-ray shell (see Hamilton & Fesen 1988); in this case, the lack of interior Ar II relative to Si II 
might be due to its much higher ionization potential (16 eV compared to 8 eV). Theoretical 
models of nucleosynthesis, accounting for heating, photoionization, and column density of 
each element would be helpful for understanding the distribution of nucleosynthetic elements. 
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The Si and S emission detected at the interior, is most hkely unshocked ejecta where 
the revere shock has not yet overtaken the ejecta. The r adial profile of unshocked ejecta 
is centrally peaked at the time of explosion, as shown by I Chevalier &: Sokerl (119891 ) . The 
radial profile of unshocke d Fe ejecta is also expe cted to be center-filled for ~1000 yr old 
Type la SNR of SN 1006 (IHamilton fc Fesenlll988l ). The morphology of the featureless dust 
resembles that of unshocked ejecta, supporting the conclusion that the featureless dust is 
also freshly formed dust. The spectrum in Figured] (curve "d"), shows the resolved two lines 
at 26 fim and at 35 fim. The two respective lines at ~26 fim may be resolved lines of [O IV] 
and [Fe II], and at ~35 /im [Si II] and [Fe II] (as expected that the unshocked ejecta near 
the explosion center have a low velocity); alternatively, they could be highly doppler-shifted 
lines (in this case the two lines at 26 fim are both [O IV], and the two lines at 35 fim are 
both [Si II]). The newly revealed unsho cked ejecta des e rves e xtensive studies; preliminary 
doppler-shifted maps were presented in iDeLaney et al.l (120061 ) and the detailed analysis of 
velocities and abundances of unshocked and shocked ejecta will be presented in future papers 
JPeLanev et allboO?! : ISmith et al l boOTh . 



3. Spectral Fitting and Dust Composition 



We performed spectral fitting to the IRS continua usin g our example reg ions in Figure 
[U Included in the fitting are MIPS 24 /im and 70 fim fluxes (IHines et al.ll2004l ). an d the con- 
tribu tion of synchrotron emission (Figs. [3] and [5]), estimat ed from the radio f luxes (IDeLaney 



20041 ) and Infrared Array Camera (IRAC) 3.6 /im fluxes ( lEnnis et al.l 120061 ). We measured 



synchrotron radiation components for each position using radio maps and assuming the 
spectral index q;=-0.71 (IRho et al.ll2003l ) where log S oc a log u. Because the full- width-half- 
maximum of 24 /im is smaller than the IRS extracted region, the surface brightnesses for 24 
fim were measured using a 15" box, the same size as the area used for the extracted IRS LL 
spectra. We also made color corrections to each MIPS 24 fim data point based on each IRS 
spectrum and band-filter shape; the correction was as high as 25% for some positions. While 
the uncertainty of calibration errors in IRAC is 3-4%, that of MIPS 24 /xm is better than 
10%. The MIPS 70 fim image (IHines et al.ll2004l ). sh own in Figure 2fe, cl early resolves Cas A 



from background emission, unlike the 160 fim image (IKrause et al. 



20051 ). Most of the bright 



70 fim emission appears at the bright ring and corresponds to the 21 fim dust map and the 
shocked ejecta, particularly [Ar II], indicating that the 70 fim emission is primarily from 
freshly formed dust in the ejecta. The 70 fim emission also appears at the interior as shown 
in Figure [2^. We measured the brightness for 70 fim within a circle of radius 20" for each 
position, accounting for the point-spread function (note that when the emission is uniform, 
the aperture size does not affect the surface brightness). We estimated the uncertainties of 
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the 70 yum fluxes to be as large as 30%. The largest uncertainty comes from background 
variation due to cirrus structures based on our selection of two background areas, 5' to the 
northwest and south of the Cas A. 

The dust continuum is fit with the Planck function {B,y(T)) multiplied by the absorp- 
tion efficiency (Qabs) for various dust compositions, varying the amplitude and temperature 
of each component. To determine the dust composition, we consider not only the grain 
species predicted by th e model of dust forrn ation in SNe (TF, N03), but also Mg protosili- 
cates (ADM) and FeO jHenning et al.lllQQsh as possible contributors to the 21 /im feature. 

culat ion are the same as those o f 
1985h. amorphous Si02 (iPhilipp 



The optical con s tants of the grain species used in the ca . 



Hirashita et al.l (120051 ) . ex cept for anaorpho us Si ( 



'iller 



19851 ). amorphous AI9O3 (IBegeniann 1 119971 ) . FeO (jHenning et al.lll995l ). and we apply Mie 
theory ( Bohren fc Huffmanf 1983 ) to calculate the absorption efficiencies, Qats, assuming the 
grains are spheres of radii a = 0.01 fim. We fit both amorphous and crystalline grains for 
each composition, but it turned out that the fit results in Cas A (see §3) favor amorphous 
over crystalline grains. Thus, default grain composition indicates amorphous, hereafter. For 
Mg protosilicate, t he absorption coe f ficient s are evaluated from the mass absorption coeffi- 
cients tabulated in iDorschner et al.l (Il980l ). and we assume that the absorption coefficient 
varies as for A > 40 /xm, typical for silicates. We fit the flux density for each spectral type 
using scale factors Ci for each grain type i, such that Fl = Sj Ci Qahs^ilO" Note that the 
calculated values of (^abs/o, are independent of the grain size as long as27r|m|a/A<<l where 
m is the complex refractive index. Thus the derived scale factor Cj as well as the estimated 
dust mass (see §4) are independent of the radius of the dust. The dust compositions of the 
best flts are summarized in Table 1. 

The strong 21 /im-peak dust is best flt by Mg proto-silicate, amorphous Si02 and FeO 
grains (with temperatures of 60-120 K) as shown in Figure [21 These provide a good match 
to the 21 /xm feature. ADM suggested that the 21 /xm feature is best flt by Mg proto-silicate 
while DOl suggested it is best flt by Si02 instead. We found, however, that Si02 produced a 
21 ixm feature that was too sharp. We also flt the observations using Mg2Si04, which exhibits 
a feature around 20 /xm and the overall variation of abso rption coefficients of Mg2Si04 with 



wavelength might be similar to that of Mg protosilicate (IDorschner et al.l Il980l : iJager et al. 



2OO3I ) . However, with Mg2Si04, the fit is not as good as that of Mg protosilicate, not only 
at the 21-/im peak, but also at shorter (10-20 /im) and longer (70 /im) wavelengths. Thus, 
we use Mg protosilicate and Si02 as silicates to fit the 21 /im-peak dust feature. The fit 
with Mg protosilicate, Si02 and FeO is improved by adding aluminum oxide (AI2O3, 83 K) 
and FeS (150 K), where AI2O3 improved the overall continuum shape between 10-70 /im 
and FeS improved the continuum between 30-40 /xm (underneath the lines of Si, S and Fe), 
as shown Figure [31 The silicate composition is responsible for the 21 /im peak, suggesting 



- 8 - 



that the dust forms around the inner-oxygen and S-Si layers and is consistent with Ar being 
one of the oxygen burning products. We also include amorphous MgSiOs (480 K) and Si02 
(300 K) to account for the emission feature around the 9.8 /xm. The composition of the low 
temperature (40-90 K) dust component necessary for reproducing 70 /im is rather unclear. 
Either AI2O3 (80 K) (Model A in Table 1) or Fe (100 K) (Model B in Table 1 and Figure HD 
can fit equally well, as listed in Table 1. We could use carbon instead of AI2O3 or Fe, but the 
line and dust compositions suggest the emission is from inner O, S-Si layers, where carbon 
dust is not expected. There are still residuals in the fit from the feature peaking at 21 /im 
(20-23 /xm), and an unknown dust feature at 11-12.5 fim (it is not a part of typical PAH 
feature), as shown in Figure [31 The former may be due to non-spherical grains or different 
sizes of grains. 

The weak 21 fim continuum is fit by FeO and Mg2Si04 or Mg protosilicate (Models 
C and D in Table 1) since the curvature of the continuum changes at 20-21 fim as shown 
in Figure [51 To fit the rest of the spectrum, we use glassy carbon dust and AI2O3 grains. 
The glassy carbon grains (220 K) can account for the smooth curvature in the continuum 
between 8-14 fim. Carbon dust (80 K) and AI2O3 (100 K) contribute to the continuum 
between 15-25 /im. We could use Fe dust instead, but we suspect carbon dust because of the 
presence of relatively strong Ne line emission with the weak 21 fim dust class. Ne, Mg, and 
Al are all carbon burning products. We cannot fit the spectrum replacing carbon by AI2O3 
with a single or two temperatures because Qabs/d of AI2O3 has a shallow bump around 27 
/im, thus the fit requires three temperature components of AI3O2 or a combination of two 
temperature components of AI3O2 and a temperature component of carbon. The continuum 
between 33-40 /im (underneath the lines of Si, S and Fe) can be optimally fit by FeS grains. 

The 70 /im image shown in Figure [2^ shows interior emission similar to the unshocked 
ejecta but that may also be due to projected circumstellar dust at the forward shock. In 
order to fit the featureless spectrum out to 70 /im, we must first correct for possible pro- 
jected circumstellar dust emission. The exterior forward shock emission is most evident in 
the northern and northwestern shell. Taking the typical brightness in the NW shell (~20 
MJy sr~^), and assuming the forward shock is a shell with 12% radial thickness, the pro- 
jected brightness is less than 4-10% of the interior emission (~40 MJy sr~^ after background 
subtraction). We assume that the remaining wide-spread interior 70 /im emission is from 
relatively cold, unshocked ejecta. Using the "corrected" 70 /im flux, the featureless spectra 
are equally reproduced by three models (Models E, F, and G) in Table 1 and Figures [61 and 
[3 All fits include MgSiOs, FeO and Si, and either aluminum oxide, Fe, or a combination of 
the two are required at long wavelength. Carbon dust can also produce featureless spectra at 
low temperature but we exclude this composition because of the lack of Ne (produced from 
carbon burning). Aluminum oxide and Fe dust are far more likely to be associated with the 
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unshocked ejecta because they result from O-burning and Si-burning, respectively and the 
unshocked ejecta exhibit Si, S, and 0+Fe line emission. However, one of the key challenges 
in SN ejecta dust is to understand featureless dust such as Fe, C, and aluminum oxide, and 
to link it to the associated nucleosynthetic products. 



4. Dust Mass 

We estimated the amount of freshly formed dust in Cas A based on our dust model fit 
to each of the representative 21 /im-peak, weak-21 /im, and featureless spectra (Fig. [1]). The 
dust mass of i-grain type is given by: 

Bu{Td,i) Ki Bu^TdA) 3 Qabs,i 

where is the flux from z-grain species, d is the distance, Bjy is the Planck function, 
is the bulk density, and a is the dust particle size. By employing the scale factor Ci and 
the dust temperature T^^i derived from the spectral fit, the total dust mass is given by 
Mdust = Yli Pi^^'^^i/'^^ where Q is the solid angle of the source. The total mass of the 21 
/im-peak dust is then determined by summing the flux of all the pixels in the 21 /im-peak 
dust region (red region in Fig. [2]F) and assuming each pixel in this region has the same dust 
composition as the spectrum in Fig. [31 We took the same steps for the weak-21 /im dust 
and the featureless dust. 

The estimated total masses for each type of dust using a distance of 3.4 kpc (Reed et al. 
1995) are listed in Table 1. Using the least massive composition in Table 1 for each of the 
three dust classes yields a total mass of 0.020 Mq (the sum of masses from Models A, D, and 
F). Using the most-massive composition for each of the three dust classes yields a total mass 
of 0.054 Mq (the sum of masses from Models B, C, and E). The primary uncertainty in the 
total dust mass between 0.020 and 0.054 Mq is due to the selection of the dust composition, 
in particular for the featureless dust. 

We also extracted a global spectrum of Cas A , but excluding most of the exterior 
forward shock regions. The spectrum is well fit with the combination of our three types 
of dust (including all compositions from Models A-G), as shown in Figure [H We used 
the dust composition of Models A-G as a guideline in fitting the global spectrum, because 
the dust features (which were noticeable in representative spectra) were smeared out. Our 
goal in fitting the global spectrum is to confirm consistency between the mass derived from 
global spectrum and that derived from representative spectra described above. The total 
estimated mass from the global spectrum fit is ~0.028 Mq, being consistent to the mass 
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determined from the individual fits to each dust class. The respective dust mass for each 
grain composition is listed in Table 3. The masses of MgSiOs, Si02, FeS and Si are more 
than a factor of ten to hundred smaller than the predictions; the predictions (N03 and TF) 
also have the dust features at 9 /im for MgSiOs, 21 /im for Si02, and 30-40 /im for FeS 
stronger than the observed spectra if the dust mass is increased. The carbon mass is also a 
factor of 10 lower than the predictions. We were not able to fit the data with as much carbon 
dust mass as expected, even if we use the maximum carbon contribution allowed from the 
spectral fits. 



5. Discussion 

We find an estimated total freshly-formed dust mass of 0.020-0.054 is required to 
produce the mid-infrared continuum up to 70 /xm. The dust mass we derive is orders of 
magnitude higher than the two previous infrared estimates of 3.5x10"'^ Mq and 7.7x 10~^ 
Mq, which are derived by extrapolation from 1.6x10""^ Mq (DOl) and 2.8x 10"^ Mq (ADM) 
for selected knots, respectively. One of the primary reasons for our higher mass estimate 
is that we include fluxes up to 70 fim while the fits in DOl and ADM accounted for dust 
emission only up to 30 and ~40 fim, respectively. The cold dust (40-150 K) has much 
more mass than the warmer (>150 K) dust. In addition, our IRS mapping over nearly the 
entire extent of Gas A with higher spatial and spectral resolutions provides more accurate 
measurements, while DOl and ADM covered only a portion of the remnant. In addition, 
ADM use only Mg protosilicate dust; the absorption coefficient for Mg protosilicate is a few 
times larger than those of other compositions. 

Our dust mass estimate is also at least one order of magnitude higher than the estimate 
of 3x 10-3 M0 bv iHines et aP tooi ). They fitted MSX and Spitzer MIPS data with Mg 



protosilicate. Note that they used only one composition. They derived a freshly synthesized 
dust mass of 3x 10^^ Mq at a temperature of 79-82 K and a smaller dust mass of 5x10^® Mq 
at a higher temperature of 226-268 K, and they explained that the mass estimate depends 
on the chosen dust temperature. As ADM mentioned, the absorption coefficient for Mg 
protosilicate is a few times larger than those of other compositions. Therefore, even including 
the long- wavelength data, the esti mated mass was sm all since only Mg protosilicate was 



modeled. With the photometry in iHines et al.l (120041 ). one could easily fit the data with 



only Mg protosilicate and would not need additional grain compositions. However, with the 
accurate IRS data, many dust features and the detailed continuum shape could not be fit 
solely with the Mg protosilicate. Note that the continuum shapes of weak 21 /xm dust and 
"featureless" dust are very different from the shape of protosilicate absorption coefficient. 
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Therefore, it was necessary to include many other compositions in order to reproduce the 
observed IRS spectra. 

It should be noted here that, in contrast with the previous works, we introduced Si and 
Fe-bearing materials such as Si, Fe, FeS and FeO. We explain why we included such dust in 
our model fitting as follows. Firstly, we inclu ded Si and Fe dust becau se these elements are 
significant outputs of nucleosynthesis; indeed IWoosley fc Weaverl (1l995l ) show that Si and Fe 
are primary products in the innermost layers of the ejecta. Secondly, we observed strong Si 
and Fe lines in the infrared and X-ray spectra; strong Si lines were detected in the Spitzer 
spectra, as shown in Figure [1] (also see DOl), and the Fe line detection at 17. 9 fim is also 



shown in Figure [5l (The Fe maps at 17.9 fim and at 1.64 /xm were presented in lEnnis et al. 



( 2006 ) and Rho et al] (12003), respecti vely.) Si and Fe lines from ejecta are also bright in X-ray 
emission (IHwang. Holt fc Petrell2000l ). Thirdly, dust such as Si, Fe, FeO and FeS is predicted 
to form in the ejecta of Population III supernovae (N03). TF and N03 predict Fe304 instead 
of FeO in the uniformly mixed ejecta where the elemental composition is oxygen-rich, but 
the kind of iron-bearing grains in oxygen-rich layers of the ejecta is still uncertain, partly 
because the surface energy of iron is very sensitive t o the concentration of impurities such as 
O and S (as was discussed by lKozasa fc Hasegawal (Il988l )). and partly because the chemical 
reactions at the condensation of Fe-bearing dust is not well understood. Depending on the 
elemental composition and the physical conditions in the ejecta, it is possible that Fe, FeO 
and/or FeS form in the oxygen-rich layers of Galactic SNe. The observations of Gas A favor 
FeO dust over Fe304, in order to match the spectral shape of the 21 /xm-peak dust and 
the weak-21 fim dust. This aspect should be explored theoretically in comparison with the 
observations in the future. 

Our total mass estimate is a lso about one order of magnitude higher than the estimate 
of 6.9x10^^ Mq by iDwekl (119871 ). who used IRAS fluxes (possibly confused by background 
cirrus) and assumed a sihcate-type dust as stellar or supernova condensates being present 
in supernova cavi ty and heated up by the reverse shock. Our estimated mass is much less 
than 1 M0, which IWilson fc Batrla I (120051 ) suggested may still be associated with the ejecta, 
after accounting for results of high-resolution GO observations. Our estimated mass of 0.020 
to 0.054 Mq is only derived for wavelengths up to 70 /xm, so it is still possible that the 
total freshly-formed dust mass in Gas A is higher than our estimate because there may be 
colder dust present. Future longer-wavelength observations with Herschel, SGUBA-2 and 
ALMA are required to determine if this is the case. Also note that we did not include any 
mass from fast moving knots projected into the same positions as the forward shock, such 
as in the northeast and southwest jets, and the eastern portions of the SNR outside the 21 
/xm-peak dust region (see Fig. [2^), because such dust could not be cleanly separated from 
the interstellar /circumstellar dust. 
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We can use our dust mass estimate in conjunction with the models of N03 and TF to 
understand the dust observed in the early universe. If the progenitor of Cas A was 15 Mq, 
our estimated dust mass (0.020-0.054 Mq) is 7-18% of the 0.3 Mq predicted by the models. 
If the progenitor mass was 30 Mq, then the dust mass is 2-5% of the 1.1 Mq predicted by the 
models. One reason our dust mass is lower than predicted by the models is that we cannot 
evaluate the mass of very cold dust residing in the remnant from the observered spectra up 
to 70 /im as described above, unless the predicted mass is overestimated. Another reason 
is that when and how much dust in the remnant is swept up by the reverse shock is highly 
dependent on the thickness of the hydrogen envelope at the time of explosion and that the 
evolution and destruction of dust grains formed in SNe strongly depend not only on their 
initial sizes but also the density of ambient interstellar medium (Nozawa et al. 2007). Dust 
formation occurs within a few hundred days after the SN explosion (Kozasa et al. 1989; 
TF; N03). Without a thick hydrogen envelope, given an age for Cas A of ~300 years, a 
significant component of dust may have already been destroyed if dust grains formed in the 
ejecta were populated by very small-sized grains; otherwise, it is possible that some grain 
types may be larger, which would increase the inferred mass. 

We observed most of the dust compositions predicted by SN Type II models, and the 
global ejecta composition is consistent with the unmixed-case N03 model than mixed-case 
model; however, note that different morphologies of Ar and Si maps imply that some degree 
of mixing has occurred. Our estimated dust mass with Spitzer data is one order of magnitude 
smaller than the predicted models of dust formation in SNe ejecta by N03 and TF, but one 
to two orders of magnitude higher than the previous estimations. We now compare the 
dust mass in high-redshift galaxies with the observed dust mass of Cas A based on the 
chemical evolution model of Morgan & Edmunds (2001). By a redshift of 4, SNe have been 
injecting dust in galaxies for over 2 billion years and there is enough dust from SNe to explain 
the lower limit on the dust masses (^7x10'^ Mf;^) inferred in submm galaxies and distant 



quasars (IChini fc Kruegel Ill994l : llsaak et al.ll2002l ). It should be noted with the dust mass 



per SN implied by our results for Cas A alone, the interpretation of dust injection from 
SNe is limited, because the amount of dust built up over time is stron gly dependent on the 



initial mass function, stellar evolution models and star formation rates (IMorgan fc Edmunds 



20031 ) . and destruction rates in supernova are believed to be important at timescales greater 
than a few billion years. Additional infrared/submm observations of other young supernova 
remnants and supernovae are crucial to measure physical processes of dust formation in 
SNe including the dust size distribution, composition and dependence on nucleosynthetic 
products and environment, and to understand the dust in the early Universe in terms of 
dust injection from SNe. 
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6. Conclusion 

1. We presented Spitzer IRS mapping covering nearly the entire extent of Cas A and 
examined if SNe are primary dust formation sites that can be used to explain the high 
quantity of dust observed in the early Universe. 

2. The ejecta maps, show a remarkable similarity to the dust maps, thereby confirming 
that dust formation occurs in the SN ejecta. 

3. The IRS spectra of Cas A show a few dust features such as an unique 21 /im peak 
in the continuum from Mg protosilicate, Si02, and FeO. We observed most of the dust 
compositions predicted by SN Type II dust models. However, the dust features in Cas A 
favour Mg protosilicate rather than Mg2Si04, and FeO rather than Fe304. The composition 
infers that the ejecta are unmixed. 

4. Our total estimated dust mass with Spitzer observations ranging from 5.5 - 70 fim is 
0.020-0.054 Mq, one order of magnitude smaller than the predicted models of dust formation 
in SNe ejecta by N03 and TF, but one or more orders of magnitude higher than the previous 
estimations. The freshly formed dust mass derived from Cas A is sufficient from SNe to 
explain the lower limit on the dust masses in high redshift galaxies. 

J. Rho thanks U. Hwang for helpful discussion of X-ray emission of Cas A. This work 
is based on observations made with the Spitzer Space Telescope, which is operated by the 
Jet Propulsion Laboratory, Cahfornia Institute of Technology, under NASA contract 1407. 
Partial support for this work was provided by NASA through an GO award issued by 
JPL/Caltech. 
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Table 1. Properties of Freshly Formed Dust in Cas A 



Dust Type 


Model 


Compositions^ 


Strong 


Nucleosynthesis 


Mass 


(spectrum in Fig. [TJ 






Lines 


Layers 


(M©) 


21/jm-pcak (a) 


A 


Mg protosilicate, MgSiOa, Si02, FeO, FeS, Si, ^^203 


Ar 


Inner-O, S-Si 


0.0030 


21/im-peak 


B 


Mg protosilicate, MgSiOs, FeO, Si02, FeO, FeS, Si, Fe 


Ar 


Inner-O, S-Si 


0.0120 


woak-21/im (b) 


C 


C-glass, FeO, AI2O3, Si, Mg2Si04 


Ne, Si, Ar (S, O+Fe) 


C-burning 


0.0180 


weak-21/jm 


D 


C-glass, FeO, AI2O3, Si, FeS, Mg protosilicate 


Ne, Si, Ar (S, O+Fe) 


C-burning 


0.0157 


Featureless (d) 


E 


MgSi03, Si, FeS, Fe, Mg2Si04 


Si, S, (O+Fe) 


0, Al burning (Fe-Si-S) 


0.0245 


Featureless 


F 


MgSiOs, Si, FeS, Fe, AI2O3 


Si, S, (O+Fe) 


0, Al burning (Fe-Si-S) 


0.0171 


Featureless 


G 


MgSiOs, Si, FeS, AI2O3, Mg2Si04 


Si, S, (O+Fe) 


0, Al burning (Fe-Si-S) 


0.0009 



^Compositions in the best fit, where a few primary compositions are written in bold, and alternative dust compositions are in italics. 



Table 2. Mass of dust with each composition 



Composition 


Model 


Model B 


Model C 


Model D 


Model E 


Model F 


Model G 


AI2O3 


6.66E-05 (083) 


O.OOE+00 (000) 


5.13E-05 (105) 


1.03E-04 (100) 


O.OOE+00 (000) 


8.13E-04 (050) 


6.50E-04 (060) 


C glass 


O.OOE+00 (000) 


O.OOE+00 (000) 


2.08E-03 (80/180) 


1.07E-03 (80/220) 


O.OOE+00 (000) 


O.OOE+00 (000) 


O.OOE+00 (000) 


MgSiOs 


1.19E-08 (480) 


1.19E-08 (480) 


O.OOE+00 (000) 


O.OOE+00 (000) 


2.55E-05 (110) 


3.19E-05 (110) 


2.55E-05 (110) 


Mg2Si04 


O.OOE+00 (000) 


O.OOE+00 (000) 


7.89E-05 (120) 


O.OOE+00 (000) 


1.72E-06 (130) 


O.OOE+00 (000) 


3.00E-06 (130) 


Mg protosilicate 


5.00E-05 (120) 


4.67E-05 (120) 


O.OOE+00 (000) 


3.77E-05 (120) 


O.OOE+00 (000) 


O.OOE+00 (000) 


O.OOE+00 (000) 


Si02 


2.23E-03 (060/300) 


1.40E-03 (065/300) 


O.OOE+00 (000) 


O.OOE+00 (000) 


O.OOE+00 (000) 


O.OOE+00 (000) 


O.OOE+00 (000) 


Si 


4.34E-04 (096) 


4.34E-04 (100) 


1.63E-03 (090) 


8.17E-03 (080) 


9.32E-04 (090) 


1.24E-04 (120) 


6.21E-05 (120) 


Fe 


O.OOE+00 (075) 


9.82E-03 (110) 


O.OOE+00 (000) 


O.OOE+00 (000) 


2.16E-02 (95/135) 


1.36E-02 (100/150) 


O.OOE+00 (000) 


FeO 


1.13E-04 (105) 


2.11E-04 (095) 


1.39E-02 (060) 


5.97E-03 (065) 


O.OOE+00 (000) 


O.OOE+00 (000) 


O.OOE+00 (000) 


FeS 


1.20E-04 (150) 


2.11E-04 (150) 


O.OOE+00 (000) 


3.40E-04 (120) 


1.94E-03 (055) 


2.59E-03 (055) 


1.29E-04 (100) 



''The numbers in parentheses are dust temperatures, and two numbers indicate two temperatures. 

I 

I— » 

00 

I 
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Table 3. Mass of dust with each composition 



Composition 


Predicted {Mq) 


Model (A+D+G)'' 


Model (B+C+E)'' 


Mass (global)'^ 


AI2O3 


2.40E-04 


~ 9.00E-03 


8.20E-04 


5.1300E-05 


1.22E-04 (105) 


carbon 


7.00E-02 


~ 3.00E-01 


1.07E-03 


2.0767E-03 


2.04E-03 (070/265) 


MgSi03 


2.00E-03 


~ 7.00E-3 


2.55E-05 


2.5500E-05 


1.65E-04 (110) 


Mg2Si04 


3.70E-02 


~ 4.40E-1 


3.00E-06 


8.0620E-05 


3.21E-05 (120) 


Mg protosilicate 


none 




8.77E-05 


4.6710E-05 


6.70E-05 (110) 


SiOa 


2.50E-02 


~ 1.400E-01 


2.23E-03 


1.3964E-03 


1.35E-03 (065) 


Si 


7.00E-02 


~ 3.00E-01 


8.66E-03 


2.9989E-03 


4.42E-03 (080) 


Fe 


2.00E-02 


~ 4.00E-02 


O.OOE+00 


3.1459E-02 


1.03E-02 (090) 


FeO 


none 




6.08E-03 


1.4136E-02 


6.23E-03 (070) 


FeS 


4.00E-02 


~ l.lOE-01 


5.90E-04 


2.1501E-03 


2.90E-03 (090) 



*The sum of the masses using the least massive composition among models. 
''The sum of the masses using the most-massive compositions. 
"^The mass using the global spectrum, 
''see the text for details. 
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Fig. 1. — Representative Set of Spitzer IRS spectra of Cas A. The contrast between contin- 
uum shape and hne emission is noticeable. Many of the strong ejecta lines are seen including 
[Ar II] (7 /im), [Ar III] (8.9 /im), [Si IV] (10.5 /xm), [Ne II] (12.8 /xm), [Ne III] (15.5 
/im), [S III] (18.7 /im), [O IV] and [Fe II] (26 /im), [S III] (33.5 /im), [Si II] (34.8 /im), and 
[Fe II] (35.3 /im). The dominant continuum shapes are 21 /xm-peak dust showing a dust 
feature at 21 fim often accompanied by a silicate emission feature at 9.8 /im with strong Ar 
lines (red, curve a) and weak-21 /xm dust with relatively strong Ne lines compared with Ar 
(green curve, b). Featureless spectra include the continuous rising spectra (yellow, c), and 
the gently rising spectra (blue, d) with strong 0+Fe or Si lines. The positions of RA. and 
Dec. are 350.900, 58.8356 (a), 350.812, 58.8075 (b), and 350.879 58.7911 (c), 350.857, 58.815 
(d), and 350.862 58.8550 9 (e). For illustration, the spectra were multiplied by 1.4 (curve 
a), 2.8 (b), 0.8 (c), 1.5 (d), and 2 (e), respectively. "Broad" continuum spectra arise from 
interstellar /circumstellar medium (orange, e). 
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Fig. 2. — (a) 21 /iin dust map; a continuum map of 19-23 /im subtracted by the baselines of 
neighboring wavelengths. This dust map is remarkably similar to the [Ar II] map (b) the res- 
olutions are convolved to match to each other). The image is centered at R.A. 23'^23™25.86** 
and Dec. +58°49'14" (J2000), and covers an 7.87' by 5' field of view, (c) A combination of 
[O IV] and [ Fe II] line map at 26 /im. (d) [Si II] (34.8 /xm) map. (e) MIPS 70 /xm map from 



Hines et al.l (120041 ). The locations of the forward and reverse shock bo undaries are marked 



as ellipses where we adj usted for elongatio n from iGotthelf et al.l (l200ll ) using long-exposure 
Chandra archival data (iHwang et al.ll2004j ). (f) Distributions of three major groups of dust 
types : 21/im-peak dust regions are in red (spectrum (a)), weak-21yum dust regions are in 
green (spectrum (b)), and featureless-dust regions are in blue ((d) spectrum in Fig. [1]). (Fig. 

2 is a jpeg file ) 
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Fig. 3. — 21 /iin-peak dust spectrum superposed on the dust fit of Model A: A Spitzer IRS 
spectrum towards a bright part of the northern shell fitted with dust compositions of Mg 
proto-silicate, MgSiOa, Si02, FeO, and AI2O3. The compositions suggest that the dust forms 
around inner-oxygen and S-Si layers. The data and the total fit are shown in blue and thick 
red lines, respectively, and MIPS fiuxes are marked with squares. The dust temperatures are 
shown in parentheses, and the dotted lines are from the second temperature components. 
Synchrotron continuum contribution (green dashed line) is estimated based on the radio 
fluxes and IRAC 3.6/xm image. 
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Fig. 4. — 21 /iin-peak dust superposed on the dust fit of Model B. 
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Fig. 5. — Weak-21 /iin dust superposed on the dust fit of Model D: A second type of dust 
continuum in Cas A. Tlie distribution of this type of dust is shown in Fig. [2]F, in green. 
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Fig. 6. — Featureless Dust Spectrum: the continuum can be fit with MgSiOs and Fe (Model 
E). The featureless spectra accompanies with S, Si, and 0/Fe lines. The green dashed line 
is predicted synchrotron emission model. 



- 26 - 




10 

Wavelength (//m) 



Fig. 7. — Featureless Dust Spectrum: The continuum can be fit with AI2O3, and MgSiOs 
(Model G). 
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Fig. 8. — Total spectrum of Cas A and the dust fit of grains with a few primary contributions 
of grain models. The dark green and brown lines are MgSiOa (100 K) and Si02 (65K), 
respectively. The green dashed line is predicted synchrotron emission model. 
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